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Abstract. The two linearly independent solutions of the three-term 
recurrence relation of the associated Askey- Wilson polynomials, found 
by Ismail and Rahman in j22|. are slightly modified so as to make it 
transparent that these functions satisfy a beautiful symmetry property. 
It essentially means that the geometric and the spectral parameters are 
interchangeable in these functions. We call the resulting functions the 
Askey- Wilson functions. Then, we show that by adding bound states 
(with arbitrary weights) at specific points outside of the continuous spec- 
trum of some instances of the Askey- Wilson difference operator, we can 
generate functions that satisfy a doubly infinite three-term recursion re- 
lation and are also eigenfunctions of g-difference operators of arbitrary 
orders. Our result provides a discrete analogue of the solutions of the 
purely differential version of the bispectral problem that were discovered 
in the pioneering work [S| of Duistermaat and Griinbaum. 



1. Introduction 

Orthogonal polynomials which are eigenfunctions of a differential oper- 
ator have a long history. When the differential operator is of order two, 
Bochner (1929) ^ proved that this property characterizes the so-called clas- 
sical orthogonal polynomials, linked with the names of Hermite, Laguerre 
and Jacobi. The general problem was raised by H.L. Krall [2HI in 1938. He 
proved that the differential operator has to be of even order and in he 
obtained a complete classification for the case of an operator of order four. 

In , it was observed that Krall's result can be reformulated in a striking 
manner. The (semi-infinite) tridiagonal matrices which encode the three- 
term recursion relation satisfied by Krall's polynomials, are obtained by 
adding a bound state (with an arbitrary weight) at the boundary of the 
continuous spectrum of the tridiagonal matrices associated with some in- 
stances of the Laguerre and the Jacobi polynomials. This way of looking 
at Krall's result leads to a very efficient construction of his polynomials, by 
means of the so-called (discrete) Darboux transformation, one of the basic 
tools in the theory of solitons. Roughly, the method consists in factorizing 
the tridiagonal matrix (appropriately shifted) as a product of a lower and 
an upper matrix, producing then a new matrix by permuting the factors. 
One of the nice outcomes of this approach is that the Darboux process can 
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be iterated, leading to orthogonal sequences of polynomials satisfying differ- 
ential equations, with a moment functional given by a weight distribution 
involving not only Dirac's delta function but also any of its derivatives. We 
refer the reader to [14| and for the precise formulation of the results in 
the context of the Laguerre and the Jacobi polynomials respectively. For 
a sample of other works related to Krall's problem or the role of the Dar- 
boux transformation in the theory of orthogonal polynomials, the reader can 
consult HZl, Cni, |2l, EZl, lini, ^ and 

This paper deals with a discrete version of Krall's problem, replacing the 
differential operator by a g-difference operator. The polynomials that one 
meets here are of more recent vintage. The celebrated Askey- Wilson poly- 
nomials [2] form the most general family of orthogonal polynomials which 
are also eigenfunctions of a second order g-difference operator, see jl2j . 
In fact, we address an extension of the problem, by allowing for a doubly 
infinite three-term recursion relation, instead of a semi-infinite one: 

The Askey- Wilson bispectral problem. To determine all doubly infi- 
nite tridiagonal matrices L, for which some family of eigenfunctions \I'(n, z) 
satisfying 

A„^'(n + l,z)+ z) + C„^'(n -l,z) = {z + z), (1.1) 

is also a family of eigenfunctions of a q-difference operator B in the spectral 
variable z, i.e. 

^ Di{z)^(n,q'z) = A{n)^{n,z). (1.2) 

finitely many i's aTL 

The motivation for allowing for doubly infinite tridiagonal matrices stems 
from the pioneering work of Duistermaat and Griinbaum [H] , where a purely 
differential version of the problem (with the tridiagonal operator L replaced 
by the Schrodinger operator) was posed and completely solved. The problem 
above was raised in by one of us in collaboration with F.A. Griinbaum, 
and solved there in the case of a g-difference operator of order two. The 
result is that, by allowing for doubly infinite tridiagonal matrices, the only 
solutions are provided by what we proposed to call the associated Askey- 
Wilson matrices. The entries of these matrices are obtained by making an 
arbitrary shift in the coefficients of the recursion relations satisfied by the 
Askey- Wilson polynomials, extending them over all integers. Precisely, 

where An, Bn,Cn, denote the coefficients of the three-term recursion re- 
lation satisfied by the Askey-Wilson polynomials. The coefficients (|1.3() 
define the so-called associated Askey-Wilson polynomials when the condi- 
tions ^'(— 1,2:) = and ^'(0,2:) = 1 are imposed in (|1.1() . The associated 
Askey-Wilson polynomials were extensively studied by Ismail and Rahman 
in the remarkable paper _^22 . However, as soon as e 7^ 0, the functions 
^(7, 2;),7 = n+e, that solve (ll.lj) and (|1.2j) . with a second-order g-difference 
operator, are not given by the associated polynomials, but rather by a two- 
dimensional space of contiguous basic hyper geometric functions, that we 
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called the Askey-Wilson functions. In these functions were constructed 
recursively in terms of an arbitrary solution of a second-order (/-difference 
equation (called there the Gauss- Askey-Wilson equation), but explicit for- 
mulas were not obtained. 

Our first result, Theorem 12.11 is to give explicit formulas for the Askey- 
Wilson functions in terms of basic hyper geometric series. For this, we use 
two linearly independent eigenfunctions of the Askey-Wilson second order 
difference operator L = La^b,c,d{l-, E) in terms of very- well-poised 8</'7 func- 
tions, that were found by Ismail and Rahman [22], The notation 
La,b,c,d{'yi E) reminds that the operator L depends on four arbitrary pa- 
rameters (that are usually denoted by a,b,c,d), 7 stands for n -|- e and E 
denotes the customary shift operator. We show that by multiplying the 
Ismail-Rahman functions by an appropriate z dependent (but 7 indepen- 
dent) factor, these functions satisfy a beautiful symmetry property, which 
we shall refer to as duality. It essentially means that the geometric and the 
spectral parameters 7 and z of the eigenfunctions are interchangeable, if we 
introduce an appropriate involution on the parameters a, b, c, d. A special 
case of this duality was considered by Koelink and Stokman 26 , in their 
study of the Askey-Wilson function transform. It is important to note that 
these authors focus on a very specific eigenfunction of the Askey-Wilson 
difference operator, which they call the Askey-Wilson function, because it 
is a meromorphic continuation of the Askey-Wilson polynomial in its de- 
gree. The same function was studied previously by Suslov jSH]) EH]) who 
established Fourier-Bessel type orthogonality relations for it. 

Sections 3 to 7 form the core of the paper. The main result can be found 
in Section 6, Theorem IH.2I and Corollary \<r>M We prove that starting with 
some instances of the associated Askey-Wilson matrices, and adding bound 
states (with an arbitrary weight) at a collection of specific points outside 
of the continuous spectrum of these operators, we can construct solutions 
of the Askey-Wilson bispectral problem, as stated at the beginning of this 
introduction. We naturally call these solutions Askey- Wilson type functions. 
When e = in we obtain in this way orthogonal polynomials which are 
eigenfunctions of (/-difference operators of arbitrary orders, thus providing 
g-analogues of Krall's orthogonal polynomials. As shown by formula (|6.18|) 
in Theorem 16.21 (see also Corollary 16. 6|) . when q ^ 1, all the added bound 
states accumulate at the boundary ±2 of the continuous spectrum [—2, 2] 
of the Askey-Wilson second order difference operator. It is interesting to 
note that most works on Krall's original problem deal with measures ob- 
tained by adding delta functions (with arbitrary masses) to the measure of 
orthogonality of some instances of the classical orthogonal polynomials, at 
the boundary of the interval of orthogonality, discarding the possibility of 
adding also derivatives of the delta function. To the best of our knowledge, 
situations involving derivatives of the delta function were first contemplated 
in ^Ij. In this case we don't have an orthogonality measure. Our result 
shows that, in the g-difference case, the bound states split up and we do get 
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an orthogonality measure for the resulting polynomials. We shall now de- 
scribe our strategy to establish Theorem 16. 21 the proof of which is prepared 
by the material exposed in Sections 3 to 5 of the paper. 

The basic technique to establish Theorem 16.21 is the so-called method 
of bispectral Darboux transformations which was developed by Bakalov, 
Horozov and Yakimov [2], and by Kasman and Rothstein 23 , in rela- 
tion with a program aiming at describing all bispectral commutative rings 
of differential operators. In the method was adapted to attack more 
systematically Krall's original problem. Though it was successful to pro- 
duce extensions of the so-called Krall-Laguerre polynomials, it led quickly 
to serious computational difficulties, when applied starting with the Jacobi 
polynomials. In the context of the Jacobi polynomials, the difficulties were 
overcome by Griinbaum and Yakimov ^3], by introducing a new idea that 
we shall explain below. Section 3 explains the method of bispectral Darboux 
transformations. Roughly, the method consists in finding the most general 
"rational" factorization of some constant coefficient polynomial in La^b,c,d 

m 

C = HLaAc,d) = lliLaAc,d - xdd) = VQ, (1.4) 
1=1 

for appropriate choices of the parameters a, b, c, d and of the bound states 
Xi to be added. The factorization involves m free parameters defining a new 
tridiagonal operator L, by exchanging the two factors 

C = HLa,b,c,d) =VQ^C = QV = f)(L). 

By a "rational" factorization, it is meant that each of the factors V and 
Q can be written as a polynomial in the Askey- Wilson difference operator 
La,b.c,d and the corresponding diagonal operator in H1.2|) . after factoring 
out some rational "function" in A^. The precise formulation is to be found 
in Theorem 13.11 The crucial Theorem 13.31 adapts the Griinbaum- Yakimov 
technology alluded to above, within the context of our paper. This theorem 
characterizes the bispectral Darboux transformations described in Theorem 
13.11 as those for which each factor V and Q is a difference operator with 
rational coefficients in q'' which is moreover invariant under the involution 

—7+1 

I(qi) = 'L^ and 1(E) = E-\ (1.5) 
abed 

Section 4 expresses some contiguous relations between the Askey- Wilson 
functions, in the language of the discrete Darboux transformation. This is 
used later in Section 6 to form in some cases the appropriate polynomial C 
in (|1.4|) that needs to be factorized, see Remark 16.41 Section 5 deals with 
the explicit computation of the kernel of the operator Q, from which the 
Darboux factorization C = VQ is performed in Section 6. It is quite re- 
markable that checking the invariance of the operators V and Q under the 
involution (|1.5|) , involves the use of classical formulas in the theory of basic 
hypergeometric series, such as the Sears' and Watson's transformation for- 
mulas. Section 7 illustrates Theorem 16. 21 on the simplest possible example, 
adding one bound state off the continuous spectrum. When e in (|1.3j) is put 
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to be zero, this example leads to orthogonal polynomials which are eigen- 
functions of a g-difference operator of order 4, providing a g-deformation of 
the original Krall-Jacobi polynomials. 

Section 8 of the paper requires some background in algebraic geometry. 
Though it played a decisive role in our research, it is not needed for the 
understanding of the earlier sections of the paper. Its aim is to provide 
some further solutions of the Askey- Wilson bispectral problem, as well as to 
take a first step in putting the problem in the more general context of the 
duality property satisfied by the Askey- Wilson functions, that we established 
in Section 2. This requires to allow in the statement of the problem itself 
for arbitrary difference as well as g-difference operators, leaving thus the 
context of orthogonal polynomials. 

In IB], we started a systematic study of bispectral commutative rings of 
difference operators, for which the dual ring is a ring of differential opera- 
tors, a problem which can be thought of as an extension of Kr all's original 
question. The techniques we used were adapted from the beautiful paper 
of G. Wilson 4r, who classified all bispectral rank 1 commutative rings of 
differential operators. By definition, the rank of a commutative ring of dif- 
ferential or difference operators is the greatest common divisor of the orders 
of all the operators in the ring. Here again the basic philosophy is "dual- 
ity" : rank 1 commutative rings of differential or difference operators enjoy 
some bispectral property, when the common eigenfunction of the operators 
in these rings possesses a symmetry which allows for an exchange of the 
geometric and the spectral variable. Ultimately, this property reflects itself 
in the geometry of the moduli space of algebraic curves which describes the 
spectrum (in the sense of algebraic geometry) of these rings. In the context 
of and JH]) the spectrum of the corresponding rings must be an affine 
irreducible rational curve with only cusp-like singularities. 

Applying the techniques of jlH] , in Proposition 18. 8| we exhibit a class 
of rational curves with double points which are the spectrum of rank 1 
bispectral commutative rings of difference operators, with dual ring a ring 
of g-difference operators. If we further impose to these rings to contain a 
tridiagonal operator, we are led to the special family of rational curves with 
double points given by 

= - 4) n - iq'-^' + Q-''^'))' n + (^''^' + 9"''^'))'' (1-6) 

with ki,l < i < g, an arbitrary collection of positive integers, which have 
been partitioned in two arbitrary sets J and K. This is the content of The- 
orem 8.1, which is the central result of Section 8. The corresponding tridi- 
agonal operators depend on g arbitrary parameters. They can be obtained 
by iteration of the Darboux transformation, starting from the operator 

E + E-\ (1.7) 

and adding bound states at the points q'^^/'^ -\- q~^^/'^ for i G J, and — (g'^*/^ + 
for i ^ K. It is well known that, when appropriately parametrized, 
the tridiagonal operators obtained in this way provide solitonic solutions to 
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the doubly infinite Toda lattice hierarchy, see HO]. Thus, we call these so- 
lutions Askey-Wilson type solitons. They are the analogues of the rational 
solutions of the Korteweg-de Vries hierarchy, which were found by Duis- 
termaat and Griinbaum 8^ to form one of the two families of solutions to 
the purely differential version of the bispectral problem. When g — > 1, the 
g-difference bispectral operator (appropriately scaled) in (|1.2j) becomes a dif- 
ferential operator. All the added bound states accumulate at the boundary 
±2 of the continuous spectrum [—2,2] of the operator (|1.7)) . The corre- 
sponding curves (|1.6jl acquire then only cusp-like singularities. We refer the 
reader to [Wl for a complete discussion of this limiting case. 

In conclusion, the solutions of the Askey-Wilson bispectral problem de- 
scribed in Theorem 18.11 as well as those described in Theorem 16.21 and 
Corollary 16.61 parallel the so-called rank 1 and rank 2 solutions found by 
Duistermaat and Griinbaum |8j, for the purely differential version of the 
problem. Whether we have found all solutions of the Askey-Wilson bispec- 
tral problem can only be pure speculation at this point. We feel that finding 
an appropriate "moduli space" of basic hypergeometric functions on which 
our result would follow from duality, looks to be an extremely interesting 
and challenging problem. Let us only mention that, even within the context 
of the Askey-Wilson polynomials, the deeper understanding of duality stems 
from affine Hecke algebraic considerations (see |33|). which were pioneered 
by Cherednik [HI, [ZI in his proof of some conjectures about the Macdonald 
polynomials. 



2. Askey-Wilson functions 



Throughout the paper we use the standard notations for basic hypergeo- 
metric series, following the book ^0] by Gasper and Rahman. In particular, 
we write 



(ai,a2, . . . ,ar;q)k = JJ(ai;g)fc, 



i=l 



with 



a; I 



and (a;g)oo = JJ(1 - ag'). 



i=0 



for products of g-shifted factorials, where < g < 1. The series expansion 



ai,a2, . 



E 

A:=0 



,ar]q)k 



{q,bi, . . .,bs;q)f, 



l+s—r 



defines the rfps basic hypergeometric series. 

We denote by E and E~^, respectively, the customary forward and back- 
ward shift operators, acting on a function = h{'y) by 

Eh.y = h-yJ^l, E ^h-y = Hy _\ . 

The Askey-Wilson difference operator La^b,c,dili E) is the second-order dif- 
ference operator 

La,b,cAl^ E) = + B^ld + CyE-\ (2.1) 
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where Id is the identity operator and the coefficients = A^-a,b,c,d, 
Ba,b,c,d, = C^;a,b,c,d are given by 

^ (1 - flbg7)(i _ acg7)(i _ Qrig7)(i _ abcdq^-^) 

^ a{l - abcdq^^-^){l - abcdq^^) ' ^ 



(l-a6cd(727-2)(i_a6cdg27-i) ' ^^'^^^ 

5-^ = a + a"^ - (yl^ + C^). (2.2c) 
Two Unearly independent solutions of the functional equation 

La,b,cA^^ E)h^{z) = {z + l/z)h^{z) (2.3) 

in terms of very- well-poised 807 series, were constructed by Ismail and Rah- 
man in They are given by 

{abq"f ,acq'^ ,adq'^ ,bcdq'^ /z;q)ao /«\^ 

r^[a,b,c,d;z) = — — ; r - x 

[bcq'^ ^bdq'^ ,cdq^ ,azq^;q)oo ^zJ 

sWj [bcd/qz;b/z,c/z,d/z,abcdq^~^ ,q~'^;q,qz/a) , (2.4) 

and 

{abcdq'^'^,bzq'^+^,czq'^+^,dzq'^+^,bcdzq'^;q)oo , 
s^{a, b, c, d; z) = — ; — - — ; — -— - — iaz) ' x 

^ {bcq'y,bdq-y,cdq'y,q'y+^,bcdzq^'r+^;q)oo 



8 



Wj [bcdzq'^'^ -jbcq"^ jbdq'' , cdq^ jq'''^^ , zq/a;q,az) , (2.5) 



where (a; ai, 02, 03, 04, 05; q, z) stands for the very- well-poised hyperge- 
ometric series 

a, qa^^"^, -qa^^'^, oi, 02, as, 04, 05 ,^ ^ 
a^/"^, —a^/'^,qa/ai,qa/a2, qa/a^, qa/a^, qa/a^ ' ' 

Let us denote by and D^^, respectively, the forward and backward g-shift 
operators, acting on a function h{z) by 

D,h{z) = h{qz) and D-^h{z) = h{z/q). 

Finally, let Ba^b,c,d{z, D^) be the Askey- Wilson second-order g-difference op- 
erator 

Ba,b,cA^^ Dz) = A{z)D, - [A{z) + A{l/z)] Id + A{l/z)D-\ (2.6) 
where 

With these notations we have 
Theorem 2.1. The functions 

Ry(a, b, c, d; z) = ,]^^°^ — '^'y(«) c, d] z) (2.8a) 
{bcd/z;q)oo 



8 L. HAINE AND P. ILIEV 

and 

S^{a, b, c, d; z) = - — ' °° 5^ (a, h, c, d\ z) (2.8b) 

[zq/a, zq/b, zq/c, zq/d; q)oo 

solve the bispectral problem 

LaAcAl, ^)^(7, z) = {z + l/z)^{^, z), (2.9) 

Ba,b,cAz, Dz)^{7, z) = A-,^(7, z), (2.10) 

where 

A^ = q-^{l-q-y){l-abcdq^-^). (2.11) 

Note that Rj{a, b, c, d; z) and Sy{a, b, c, d; z) differ from r^(a, 6, c, d; z) and 
s-^(a, b, c, d; z), respectively, by a factor independent of 7, hence p.9j) is auto- 
maticahy satisfied. The proof of ()2.1U|) is based on an interesting bispectral 
involution. As we shall see, when appropriately normalized, r^{a,b,c,d; z) 
and s^(a, b, c, d; z) satisfy an important duality relation. All this is the con- 
tent of the next lemma. 

Let us introduce "dual" parameters and variables via the formulas 



^ 6=4 a=#^. <i=4. (2.12a) 

abed aab aac aad 

,= C T^'f^-l, (2.12b) 



and define 



.(^ h r d-7\- /,p^,-.2-los(abcd)/ log(q) {bc, bd, cd, qz/a, qz/b, qz/c, qz/d; q)^ 

[bz, cz, az, bed/ z] q)oo 

(2.13) 

/ , , V bcdz 
.-(a,b,c,d;z)= . (2.14) 

Denote by R and S the functions 

Rj{a,b,c,d; z) = f{a,b,e, d; z)r^{a,b,e, d; z), (2.15) 
S^{a, b, c, d; z) = s(a, 6, c, d; z)s^{a, b, c, d; z). (2.16) 

Lemma 2.2. The functions R and S defined above satisfy the duality rela- 
tions 

R^{a, b, c, d; z) = R;y{d, b, c, d; z), (2-17) 

and 

Sj{a,b,c,d; z) = Sj{d,b,c,d; z). (2-18) 



Proof of Lemma \2JA. Both (|2.17j) and (|2.18)) can be proved by using trans- 
formation formula jlOj (III. 24), p. 243]. Indeed, replacing a, 6, c, d, e, / by 
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bcd/qz, abcdq"^^^, c/z, d/z, b/z, q"'^ we have 
^Wj {bcd/qz; b/z, c/z, d/z, abcdq^~^ , q~^; q, qz/a) 

{bed/ z, abcdq"^^^ , qb/a, qc/a, qd/a, zq~^~^^ /o.;q) 



— X (2.19) 

{be, bd, cd, bcdq'y /z, q~'^+'^ /a?, qz/a; q)oo 

%Wj {q^^^^ /a^; q^^^^ /ab, q^^^"^ /ac, q~^^^ /ad, qz/a, q/az; q, abcdq^^^) . 

The series in the right-hand side of (|2.19p can be rewritten in terms of 
the dual parameters (see ()2.12|) ) as 

%W^ (bcd/qz;b/z,c/z,d/z,dbcdq^~^,q~'^;q,qz/a^ , 

which is exactly the factor in r^(a, 6, c, d; z). Now taking into account 
(fTl)) . (Einj), ((TT^ and (HHH), one can easily check that 

r^{a,b,c,d; z) f{a,b,c,d;z) 
r;y{d,b,c,d;z) r{a,b,c,d, z)' 

which combined with H2.15() gives H2.17() . The proof of H2.18() is similar. 
First we apply dni (III. 24), p. 243] with a, b, c, d, e, f replaced by bcdzq^^ , 
q'^~^^, bcq^ , bdq"', cdq^ , zq/a to obtain 



^W'i {bcdzq^'^; bcq^ , bdq^ , cdq'^ ,q^~^^ , zq/a; q, az) 

{bcdzq^'^'^^ , q^~^^ , abczq'^ , abdzq^ , acdzq'^ , qz"^ ; q) 
{bzq^~^^, czq^^^ , dzq'^^^ , abcdq^'^ , abcdz"^ q^ , az; q 
8W7 {abcdz'^q'^"^; az, bz, cz, dz, abcdq"^'^ ; q, q'^'^^) ■ 



X (2.20) 



Again the hyper geometric series in the right-hand side can be written in the 
dual parameters as 



sWj {icdzq^^'; bcq^ , bdq^ , cdq^ , q'^''^^, zq/d; q, dzj , 

and the proof follows by combining ((731), 111201), 11112), (tTTll) and (171^ . 

□ 

Proof of Theorem \2.1\ Let F = F^{a, b, c, d; z) be R-y{a, b, c, d) or S-y{a, b, c, d). 
From Lemma 12.21 it follows that F satisfies also a difference equation in 7, 
i.e. we have 

^a,6,c,d(7,^7)^7(«'^'C,d;z) = {z + l/~z)F^{d,b,c,d;i). (2.21) 

On the other hand, from (|2.12bl) . it is clear that 

F^±i{d,b,c,d;z) = Df^F^{a,b,c,d; z). (2.22) 

Combining (|TT|) . and (fT^ we see that F^{a, b, c, d; z) satisfies 

the following g-difference equation in z 

d{z + 1/z — a — \/a)Fry{a, b, c, d; z) 

= A{z) {F^{a, b, c, d; zq) - F^{a, b, c, d; z)) (2.23) 

C{z) {F^{a, b, c, d; z/q) — F^{a, b, c, d; z)) , 
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where A{z) = dA!y{d, b, c, d) and C{z) = aC;y(a, b, c, d). Using ()2.2a|) . (|2.2b|) . 
and (|2.12|) . one can express the coefficients of the g-difference operator and 
the eigenvalue in formula (|2.23j) in terms of the parameters a, b, c, d and of 
the variables z and 7 as 

7/ N _ (1 - qz/a)il - qz/b){l - qz/c){l - qz/d) 

^(^) - il-z^){l-qz^) ' ^^-^^^^ 

ru \- 9^ (1 - az/q){l - bz/q){l - cz/q){l - dz/q) 

and 

d{z + l/z-d- 1/d) = q^+^ + ^ - -f!^ - 1. (2.25) 

abed abed 

If we define $(z) to be the function 

^-2+\oR{abcd)/\oR{Q) {az,bz,cz,dz;q)oo 

{zq/a, zq/b, zq/ c, zq/ d; g)oo ' 

then from (IT^ . (IT^ and TTM it follows that G{z) := ^{z)F^{a, b, c, d; z) 
satisfies the g-difference equation 

+ ^ - - 1) G{z) = A{z)G{zq) - B{z)G{z) + C{z)G{z/q), 
abed abed j 

where 









(2.27) 


A{z) 




q^ 

abed ' 


(2.28a) 


C{z) 




= -j-Aii/z), 

abed 


(2.28b) 


B{z) 


= A{z) + C{z), 




(2.28c) 



with A{z) defined in (|2.7j) . A direct computation now shows that 

A/ N ^ ( A, ^ A,^ I \ . fl&cd abed\ /„ „„x 

B{z) = 4— A[z) + A{\lz) + g - 1 + , (2.29) 



abed \ q 

which combined with (|2.27|) and (|2.28p proves that G{z) satisfies the Askey- 
Wilson q-difference equation 1)2. 10() . To finish the proof, it is enough to 
notice that 

(^{z)Ry{a, b, e, d) = bcd{bc, bd, ed; q)oo R-yia, b, c, d; z), (2.30) 

and 

^{z)S-y{a,b,e,d) = bed S^{a,b,e, d; z), (2-31) 

i.e. ^{z)R and ^{z)S differ from R and S, respectively, by unessential 
constant factors. □ 
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3. BiSPECTRAL DARBOUX TRANSFORMATIONS 

We denote by 

B={La,b,c,d.^l). (3-1) 

the algebra of difference operators generated by La^b,c,d and A^, defined by 
(|2.1j) and (|2.11j) . respectively. Similarly 

B' = {z + z-\Ba,b,c,d) , (3.2) 

will denote the algebra of g-difference operators generated by the operator of 
multiplication by z + z~^ and the operator Ba^b,c,d defined in (|2.6|) . Formulas 
(|2.9j) and (|2.1()j) serve to define an anti-isomorphism 

b-.B^B' (3.3) 

between these two algebras, i.e. it is given on the generators by 

b{La,b,c,d) = z + z^'^ and b(A^) = Ba^b,c,d- (3.4) 

Note that the bispectral property (|2.9p . l|2.1Ujl is equivalent to the identity 

x^{j, z) = b{x)^{j, z), yx e B. (3.5) 

We shall also need the commutative subalgebras (the algebras of "func- 
tions" ) of ;S and B' defined by 

IC = {Ay) and JC' = {z + z~'^). (3.6) 

The next theorem summarizes the technology of bispectral Darboux trans- 
formations, initiated in i^, [Ij and |2S|- It was adapted and applied to the 
case of difference operators in |14| . For the convenience of the reader, we 
include the short proof of this result. 

Theorem 3.1. Let C be a constant coefficient polynomial in La^b,c,d; which 
factorizes rationally as 

C = rQ, (3.7) 

in such a way that 

r = ur-^, Q = e-W, (3.8) 

with U,V € B, and 0,r € /C. Then the Darboux transform of C given by 

C = QV (3.9) 

is a bispectral operator. More precisely, defining f = h^C) G K,' and ^ = 
Q^, with ^ satisfying ((T^ and (HHUl), we have 

Ci> = f^, (3.10) 

M = er^, (3.11) 

with 

B = b{V)b{U)f~^. (3.12) 
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Proof. From the definitions above and equation (jSISJ), we obtain 

£i> = QVQm = QC^ = Qb{C)^ = fQ^ = fii, 
which estabhshes (|3.1flj) . Using (|3.5|) and (|3.8j) . we can write ^ as 

4> = Q^ = e^V^ = e"^b(y)^. (3.13) 

From 1)3 .Tp and H3.8|l . we deduce that 

er = vc-^u. 

Applying the anti- isomorphism b to this equation, we obtain 

b(er) = b{U)f'^b{V). (3.14) 
Finally, using ((231), (PO^ . ((mi) and (1X111) we get 

er^- = GrG-^b(y)^ = e"^b(y)Gr^' = e~^b(y)b(er)* 

= e-^b{V)b{U)f-^b{V)'^ = b{V)b{U)r^e-^b{V)'^ = B^, 
which gives (|3.11() and completes the proof. □ 

Despite the apparent simplicity of Theorem 13. 11 it is a priori very compli- 
cated to recognize that an operator admits a rational Darboux factorization 
as defined by ((3.7() , (|3.8p . The rest of the section is concerned with making 
Theorem 13.11 more effective. Some of the ideas developed below were first 
introduced in 15;, in the context of the Jacobi polynomials. 

Let us denote by TZ{q'~' ,E} the algebra of difference operators of the form 

7712 



j=m.i 

with coefficients hj{q^) rational functions in . The ordered pair [mi, 7712] 
is called the support of T. Consider the involution^ I on lZ{q"' , E}, defined 
by 

—7+1 

nqi) = 'L^ and 1(E) = E'K (3.15) 
abed 

Lemma 3.2. A Laurent polynomial in q"' 

n 

p{q')= "^."^N, (3.16) 

k=—m 

is a polynomial in (i.e. it belongs to K, as defined in (|3.6() ) if and only if 
it is I-invariant. 

Proof. A straightforward computation from (|2.11() gives 

/(A^) = A^, 



^By involution we mean an automorphism /, such that / o / = Id. 
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hence every polynomial in is /-invariant. Conversely, p{q"') in (|3.1(i|) is 
/-invariant, one checks that m = n and 

Pm = ^4^'^^f^)- (3.17) 



A:=0 

Observe from 1)2. 11(1 that 



-Y (7 / , ^ abed 



abed abed V q 
Since 



{abcd)^ V abed J \ {abed) 



,(fc-2)7 



,(fc-2)(l-7) 



+ 



o6cd V {abed)^ ^ 

by induction, it follows that p{q^) in ()3.17|) is a polynomial in A^, i.e. it 
belongs to /C. □ 

An easy computation, using formulas (|2.2a|) and 1)2. 2b|) . gives 

I{A.y-a,b,c,d) = C^-^a,b,c,d, (3.18) 

which shows that the operator La^b,c,d ™- (|2.1|) is /-invariant. Let us denote 
by the subalgebra of TZ{q"' ,E^, consisting of /-invariant operators, i.e. 

= {T G 7^{g^, E]\ I{T) = T] . 

Since La^h,c,d G and A^ G A^, it follows that 

ICcBcA^. (3.19) 

Let A^ be the subset of A^, consisting of /-invariant operators with support 
[—m,m]. The next theorem characterizes the bispectral Darboux transfor- 
mations as those for which the two factors V and Q in ()3.7|1 are /-invariant 
operators, with rational coefficients. 

Theorem 3.3. The following conditions on an operator T G TZ{q''',E} are 
equivalent: 

(i) The operator T is I -invariant, i.e. T G A^; 

(ii) T has the form Q~^V, for some operator V € B and some function 
e G /C; 

(iii) T has the form UT^'^, for some operator U ^ B and some function 

r G /c. 

Proof. From (|3.19|) it is obvious that (ii) implies (i). Below, we shall prove 
that (i) implies (ii), i.e. if T G A^, then T = Q~^V , for some G /C and 
V €z B. The proof of the fact that (i) is equivalent to (iii) is similar. 

It is clear that T must have support [— m, m], for some m G N. If m = 0, 
we can always assume that T is the quotient of two /-invariant Laurent 
polynomials in q"'. Hence, by Lemma l3.2| T is a rational function in A^, 
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SO we can take V and Q in (ii) to be the numerator and the denominator, 
respectively. If m > 1, we shall show that 

T = e~V„ modA^^.i, (3.20) 

with @m G ^1 Vm G S, from which the proof follows by induction. 
Any operator in may be written 

r=(Id + /)^i?- modA^;„_i, (3.21) 

with c{q'^) and (}{q''') polynomials in q"^ . In particular, since La^b,c,d is I- 
invariant, we have 

LT,,,,,,= {ld + I)^^E"^ modA^_i, (3.22) 

for some polynomials fiq"^) and siq'^)- 

We denote by A^; and px the operators of left and right multiplication by 
X, respectively, on Tl{q'^,E}, i.e. for T G TZlq^^^E} we have 

Xx{T) = xT and px{T) =Tx. 

Eliminating q~'^ (resp. q^) from the equations 

K^E"^ = (a6cV~^ +q-"' -I- abcdq-^)E"', 

E'^K^ = {q'^ahcdq'-^ + q-'^q'^' - 1 - abcdq-^)E"' , 

one obtains 

S+^E"") =q^E"' and 6-{E'^) = q-^E"', (3.23) 



with 



and 



•5^^ = . :-2n^ (V - '?">A,) + f ^ + 1 ) Id 



-m 

abed 



^ S^m- (3.24b) 

Notice that 6^ and 5~ preserve the subalgebra A^. Hence, from ()3.22|) . 
using the properties 1)3. 23() of 6^ and (5~ , one deduces that 

p(<5+ , 5- )L-,,,,, = (Id + I)piq\ Q-'')^E"^ mod A^,,, 

for any polynomial p(g''', q^'^) G c/^^]. In particular, if we choose 

Piq^q-'') = ciq-y) Q{q^) I{d{q^) f{q^)), 

we get 

p{6+,5-)L^^,^,^, = Hq^) m I{l>iq^) (Id + I)^i?- mod A^.,. 
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Define 

Since Qm is /-invariant, by Lemma 13.21 @m £ fC. Also, it is clear from 
the definitions H3.24a|) and ()3.24b|) of (5+ and 5;; that p{5m^^m)I^'^,b,c,d ^ ^■ 
Hence, remembering (|3.21|) . we get 



pi^L ^;n)La,b,c4 = (m + i)'^^e^ mod 

= T mod 

which establishes (|3.20|) with Vm = p{^m->^m)^^b cd- This completes the 
proof of Theorem 13.31 □ 

4. Contiguous relations 

Prom ()2.5() . we observe that the function s^{aq,b,c,d;a) has a simple 
product form, since the series defining it reduces to 1. One checks 
easily that 

s.y{aq,b,c, d;a) _ ip^ C^;a,b,c,d 
s^-i{aq,b,c,d;a) ^P-y-i A^.a,b,c,d' 

with 

f-y = \ , • (4.1) 

q-y - abcdq'y ^ ' 

Let us denote by Pa,b,c,d and Qa,b,c,d the difference operators 

Pa,b,c,d = ^^iE-ld), (4.2) 
Qa,b,c,d = {^■y;a,b,c,d^d. — C^.a,b,c,dE . (4.3) 

From ((TT|) and ((T^ . it follows that 

La,b,c,d - {a + a^^) Id = Qafi,c,dPafi,c,d- (4.4) 

Since the kernel of the operator Qa,b,c,d is generated by s^{aq,b,c,d;a), it 
is natural to expect that interchanging the factors Qa,b,c,d and Pa,b,c,d in 
1)4. 4(1 . we obtain the operator Laq^b,c,d- This simply means that La^b,c,d can 
be obtained as a Darboux transformation from the operator Laq^b,c,d- The 
proposition below makes this statement precise. 

Proposition 4.1. With the notations above, we have 

Laq,b,c,d — {a + a "'^) Id = Pa,b,c,dQa,b,c,d- (4.5) 

Moreover, up to a factor independent of Qa,b.c,d maps r^{aq,b,c,d; z) 
and s^{aq,b,c,d; z) into r^{a,b,c,d; z) and s^{a,b,c,d; z), respectively. More 
precisely, the following contiguous relations 

Qa,b,c,d r^{aq, b, c, d; z) = — - — r^(a, b, c, d; z), (4.6a) 

Qa,b,c,d s-yiaq, b, c, d; z) = — — '^^^ — ^s^(a, b, c, d; z), (4.6b) 

hold. In particular, from ()4.6b|) . it follows that the kernel of Qa,b,c,d is 
spanned by s^{aq,b,c, d; a) . 
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Proof. A straightforward computation shows that 

a + a~^ —C^+l-a,b,c,d — A^;a,b,c,d = aq+{aq) — A^.aq,b,c,d — C^;aq,b,c,d, (4.7) 

which combined with ((TT|) . (EHJ, and (|0)) gives The proof of 

1)4. 6() can be easily extracted from the contiguous relations found by Ismail 
and Rahman in Let us denote for simplicity by x^{a,b,c,d; z) and 

5^{a,b,c,d; z) the factor on the right-hand side of formulas (|2.4() and 
1)2. 5() . respectively. Applying formula (2.3), p. 207], with a, b, c, d, e, / 
replaced by bcd/qz, abcdq"^, q~^ , b/ z, c/ z, d/z, we obtain 

{q~'^ — abcdq^^^) ^1 ^—^ ^-/{a, b, c, d; z) 

= {q-^ - 1) (l - r,_i (ag, 6, c, d; z) (4.8) 

+ (l — abcdq^^^^ ^1 t^(ag, 6, c, d; z). 

From ()2.4|) and (|4.8() , one can deduce ()4.6a() . To get 1)4. 6b|) , we first substitute 
6cdzg27+i, z/a, bcq^ , bdq^, cdq"/ , q^+^ for , A/\ Aj [i, A/v, A/ p, A/a in 
m formula (2.17), p. 210] to obtain 

{l-bcdzq^^){l-bcdzq^^-^) 

(1 — z/a)(l — az)5-y_i(a, 6, c, d; z) (4-9) 

z [l - abq'^){l - acq^){l - adq^){l - abcdq^-^) 
+ a (l-a6cV7-i)(l-a6cdg27) ^' 

which combined with (|2.5|1 gives (|4.6b|) . □ 

Remark 4.2. We can also consider Laq^b,c,d as obtained by means of a 
Darboux transformation from La^b,c,d- In this case, the analogous formulas 
to (jTBal) and ()4.6b|) are 

Qj — Z 

Pa,b,c,d (a, b, c, d; z) = ^ {aq,b,c,d; z), (4. 10a) 

-Pa,fe,c,d ^7(0, c, d; z) = -s^(ag, 6, c, d; z). (4.10b) 
One way to prove these formulas is to use [HI formula (2.2), p. 207] and 



formula (2.18), p. 210], and to proceed as before, adapting the proof of 
Proposition 14. II However, they can be easily deduced from Proposition 14. 11 
by applying Pa,b,c,d to both sides of (|4.6a|) - ()4.6b|) . and by using (|4.5jl and 
the fact that r^(a, b, c, d; z) and s^(a, 5, c, d; z) satisfy (|2.9j) . 

5. Computing the kernel of the Darboux transform using 
Sears' and Watson's transformation formulas 

As we saw in Theorem 13.31 in order to apply the general bispectral Dar- 
boux technology, one needs to find an /-invariant factorization of a constant 
coefficient polynomial Cm Lab cd- In this section we construct /-invariant 
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eigenfunctions of La^b,c,d, which generate the kernel of the operator Q, in the 
notations of Theorem 18 .11 



Proposition 5.1. Assume that 

a = dq"' and Zk = dq'^, fc, a G N, < A; < a - 1. 



(5.1) 



Then, 



r^{zk) = {bcq ^ ,d^q°';q)k^ 



'd^q^,hcd^q^^'^~^,q~'^,q~^ 

bd, cd, d^q"" ' ^' ^ 



s^izk) = {hcr-''-\dWdq^-'^ /h,dq^^^ /c-q)^.k-l^ 



^k-a+1^ 6cgT, g-T-"+V(i2, q^-»-''/d'^ _ 



'/d\bq^-'^/d, cq^-'^/d 



with 



ibdq'y,cdq^;q)a 



g-y = q 



{q^+\bcq^;q)oo 



If we make the further assumption that 

d^ = 

g-y in H5.5() can be rewritten as 



d^ = q'', len, l> 1, 



■yil+a-l) 



9-1 



{q'f+\bcq'y;q)i+a-i' 



(5.2) 
(5.3) 

(5.4) 
(5.5) 

(5.6) 
(5.7) 



Proof. Applying Watson's transformation formula |1UI (III. 18), p. 242] with 
n, a, b, c, d, e replaced by k, bcq~^~^ , bq'^'/d, cq~^/d, bcd"^ q'^^"^^ , q^'^, 
respectively, we can rewrite r^[zk) as 



^ (,_,H (rf^g"+^feg-^g)fc (g-"-^+Vd^;g), 
^ {bdq\cdq^-q)c. {q-^—'^+^/d^;q)k 

'd^q\bcd'^q'^+''-\q-^,q-'' _ 



103 



■,q,q 



{51 



cd, bd, d'^q" 

Using formula (1.9) on page 233 in one can easily deduce that 
(g-"-fe+Vd2;g)fc ^,JdVl^ 

which shows that (|5.8j) is equivalent to (|5.2|) . The proof of l|5.3|) is similar. 
Let us first replace the parameters n, a, b, c, d, e in Watson's formula lilOl 
(III.18), p. 242] by a - A; - 1, bcd^q^^+^, bdq^ , cdq^ , hcq^ , q^+^, to obtain 
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the following formula for s^i{zk) in terms of 4(^3 series 



493 



q^^^ ,bcq^ ,q"<^^ ,q~°'^'^'^^ 
cdq^+^+\bdq^+^+\q^-^/<f''^'' 



(5.9) 



Now replacing n, a, 6, c, d, /bya— /c— 1, 6cg''', g'^"'""^, q"'~^^ , q^ '^/d'^^ cdq"'^^'^^, 
bdq^~^^~^^ in Sears' transformation formula (see |in| (III. 15), p. 242]), we see 
that ((SHI) gives exactly □ 

In the case considered above (i.e. a = dg" and d'^ = q^) formula (|3.15|) 
reduces to 

-7-a+l „-7-Q-«+l 

Lemma 5.2. If (|5.6() holds, the functions r^{zk) and s^{zk), definedby 1)5. 2() 
anrf (|5.3() are I -invariant rational functions in q^ , i.e. we have I{r^{zk)) = 
T^{zk) and I{s^{zk)) = s^{zk). 

Proof. From (|5.2|1 . 1)5. 3(1 . ()5.4p . and (|5.7() . it is obvious that r-y{zk) and 5^(2;^) 
are rational functions in q"*. Since I{bcd'^q'^~^'^~^) = q^"^ and I(bcq'~') = 
q-'^-a+i I ^ ^ it is clear that the 4(^)3 series in (|5.2() and 1)5. 3|) are /-invariant. 
Thus, to show the invariance of r^{zk) and s-y(^fc), it is enough to show the 
invariance of /-y and g-^. Using the definitions (|5.4|1 of f^ and (|5.in|l of /, 
one computes 

A/7)- ^ ^ 



(1 - l/cdg')'+°-i)(l - l/c(igT+«-2) • • • (1 - l/cdq^) 

1 

^ (1 - l/6(igT+a-i)(i _ l/6dg7+"-2) . . . (1 _ l/6dg7) 

„27Ci+ci{o—l)— 0(7+0— 1) 
= ^ = f 

{bdq^,cdq^;q)a 

The proof of the /-invariance of (7^ is similar. □ 



6. AsKEY- Wilson type functions 

In this section we prove the main result of the paper. First we review 
some basic facts about Darboux transformations. 

6.1. Iterating the Darboux transformation. Consider a difference op- 
erator Lq and a nonzero eigenfunction V'i(7) foi' which LoV'i(7) = a^iV'i(7)- 
If Qq is a first-order difference operator with kernel spanned by ^'1(7); we 
can write Lg in the form Lq = xiId + PoQo) fo'^ some difference operator Pq- 
The operator Li = xild + QqPq is by definition a Darboux transformation 
of Lq. 
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The next proposition, which can be found in [4:'A\ pp. 14-19] for the case 
of differential operators, describes the resuh of m successive Darboux trans- 
formations starting from Lq 

Lq = Xjjid + PqQo r\ Li = Xj^^Id + QqPo = Xj^ld + PiQi r\ ■ ■ ■ 

Id + Qm-2Pm-2 — + Pm-lQm-1 (6.1) 

r\ L = Lm = Xj^ld + Qm-lPm-l- 

Proposition 6.1. // the operator L is obtained from Lq by iteration of the 
Darboux transformation (|6.1j) then 

LQ = QLo, (6.2) 

where 

Q = Qm-iQm-2- ■ -Qo- (6.3) 
Conversely, if L and Lq are difference operators and if there exists a 
difference operator Q of order > 1 such that (|6.2)) holds, then L can be 
obtained from Lq by a sequence of Darboux transformations. 

Proof. From (|6.1j) one can deduce that 

LkQk~i = Qk-iLk~i for = 1, 2, . . . , m, (6.4) 

which easily implies ()6.2|) . thus proving the first assertion. 

The second part can be proved by induction on the order of the operator 
Q as follows. Assume that ()6.2() holds. It is obvious that Ker Q is preserved 
by Lq, i.e. Lo(Ker Q) C Ker Q. Thus we may regard Lq restricted to Ker Q 
as a linear operator acting in a finite-dimensional complex vector space. Let 
{V'i(7)}j=i be a Jordan basis of this operator, i.e. 

KerQ = Span{V'i(7), V'2(7), • • • , V'm(7)}, (6-5a) 

where 

^oV'fc(7) = Xj^Tpkil) + o"fcV'fc-i(7), I <k <m, (6.5b) 
with 0"! = and Ufc = or 1 for 2 < k < m. Let Qq be a first-order difference 
operator with kernel spanned by ■0i(7). If the order of Q is 1, we simply 
take Qq = Q. Otherwise we can take for example 

Since V'i(7) is in the kernels of Lq — xj-^ld and Q, there exist difference 
operators Pq and Q such that 

Lq - xjjd = PqQo (6.6) 

and 

Q = QQo. (6.7) 
If we define Li = Xj^ld + QqPq, then from (|6.6() we get that 

LiQo = QqLq. (6.8) 
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Using the last equation, (|(i.2|) and (|(i.7j) one can show that 

LQ = QLi. 



Moreover from ()6.5b|) and (|6.8|) we can see that for A; = 2,3,...,mwe have 

ii(QoV'fc) = Xj^Qoipk + a'kQoipk-i, where a'2 = 0, and = cjfc if A; > 2, 
i.e. {QoV'i(7)}jL2 is ^ Jordan basis for the operator Li in the vector space 



Ker Q. The proof now follows by induction. 



□ 



The above proposition tells us that the Darboux process (|6.H) is deter- 
mined by the intertwining operator Q in (|6.2|) . The points {a^j^} are simply 
the eigenvalues of the operator Lq in the vector space Ker Q. We will use 
Proposition 16.11 as follows. First, we will construct a difference operator Q 
whose kernel is Lq invariant, with Lq = La^b,c,d the Askey- Wilson operator 
defined in Section |2 This operator Q will correspond to a sequence of Dar- 
boux transformations (|6.1() . Indeed if we consider the operator QLq, then 
its kernel contains the kernel of Q, so it has a factorization of the form 

LQ = QLq, 

for some (unique) difference operator L. By the above considerations, L is 
obtained from Lq by a Darboux chain (|6.1j) . 

The next important fact is that the operator Q needed above can be 
explicitly reconstructed (up to a factor on the right) using the functions 
{tpj} from its kernel (|6.5a|) . More precisely, if [mi, 7712] is the support of Q, 
with rn-2 — rrii = m we can write the operator Q in the form 



Q = a(7) 



(6.9) 



-0i(7 + mi) ••• ilJm{l + mi) 

where 0(7) is some function in 7, and the determinant is expanded from 
left to right, i.e. the shift operators E^ are pulled to the right. Indeed, 
it is obvious that the right-hand side of ()6.9() defines a difference operator 
with support [mi,m2], whose kernel coincides with the kernel of Q. But 
two operators having the same supports and kernels must differ by a factor, 
which proves ()6.9|) . 

Another useful way to describe the Darboux process ()6.1() is to define a 
new operator C by the following formula 



£ = f)(Lo) = n(Lo-x,Jd). 



Using H6.1|) - ()6.4() one can easily show by induction that 

C = VQ, 

where 



V = PqP,--- r 



m—1 ■ 



(6.10) 



(6.11) 



(6.12) 



ASKEY-WILSON TYPE FUNCTIONS, WITH BOUND STATES 21 
Thus, from (|(i.2j) we get 

m 

C:=QV = J{{L-x,ld) = \){L), (6.13) 

i=l 

i.e. we have the Darboux map 

C = \){Lq)=VQ^C = QV = \){L). (6.14) 

6.2. The main result. Now we are ready to formulate our mam result. 
We show in Theorem 16.21 below, that if 

d = ibg'/^ and a = dq", for some Z,a G N, l,a >1, (6.15) 

then, the second-order difference operator L obtained by iteration of the Dar- 
boux process ()6.1|) . starting with the Askey- Wilson operator Lq = La^b,c,d, 
and iterating at any subset {xj-^ , ■ ■ ■ , Xj^} of the points 

Xfe = zfc-Fz^\ Zk = dq^ = ±q^/'^'^^ 0<k<a-l, (6.16) 

possesses a two-dimensional space of eigenfunctions which are also eigen- 
functions of a g-difference operator in the spectral variable z. 

Notice that in this case the functions {ipj} must be linear combinations 
of the functions rj{zj^^) and s-y{zj^), computed in Proposition 15.11 i.e. 

i^kil) = tJ'kr'y{a,b,c,d;zj^) + h'kS^{a,h,c,d\ Zj^), k = 1,2, . . . ,m, (6.17) 

for some constants {fJ'kT^k}k=i- Moreover the operator Lq is diagonalizable 
(in the vector space KerQ), i.e. = for all k in (|6.5b|) . It is interesting 
that this is a purely g-phenomenon. In the q = 1 case, one needs to use 
generalized eigenvectors, see ^S] and flU] . 

Theorem 6.2. Assume that ()6.15|) holds, i.e. a = ±52+° and d = ^q^ , 
l,a €z N, l,a > 1, with the same choice of sign for a and d. The second- 
order difference operator L which is obtained by iteration of the Darboux 
transformation ()6.1() . starting with the associated Askey-Wilson operator 
Lq = Lafi^cd (|2.1|) . and adding bound states at the points 

{xfc = dq^ + {dq^)-^ = ± [q^+^ + g-(5+^0) , 

k G {ii, j2, . . . , i„^} C {0, 1, . . . , a - 1}}, (6.18) 

possesses a two-dimensional space of eigenfunctions which are also eigen- 
functions of a q-difference operator in the spectral variable z. 

If ^'(7,2:) satisfies (|2.9() (i.e. it is a linear combination of and 
defined by (|2.8jn . then the function ^(7, z) defined by 

i>{-i,z) = Q^>{-i,z), (6.19) 

satisfies 

L|f(7, z) = {z + z-^)^{-i, z). (6.20) 
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Equation (|().2())) follows immediately from (|().2|) . So, it remains to show that 
there exists a g-difference operator B[z,Dz) in z, such that 

^(z,Z),)l'(7,z) = A^4'(7,z), (6.21) 

for some function A^. 

Before we start with the proof of this fact, we shall make two simple but 
important remarks. 

Remark 6.3. So far, we have not said anything about the normalization 
and the supports of the operators {Pk,Qk}^=Q ^ or equivalently, about the 
normalization 0(7) and the support [mi,m2] of Q in (|6.9() . The point is 
that this is not important for the bispectral relation 1)6.21(1 . Indeed, if H6.21() 
holds, and if 

q! = a'(7)^™' Q and l''(7, z) = Q'^{-/, z) = a'(7)l'(7 + m', z), (6.22) 

for some function o'(7) and some m' G Z, then it follows from (|6.21() and 
((?r^ that 

B{z,D,)i>'{-i,z)=k'^^>'{^,z), (6.23) 

where k'^ = A^+m', i.e. we still have a g-difference equation in z. Thus, the 
choice of the support [mi, 7x12] and of the normalizing function 0(7) in (|6.9() 
is not essential for the bispectrality. 

Remark 6.4. From Proposition 14.11 we know that -La,fe,c,d is a Darboux 
transform from Laq^b,c,d at a+a~^. Thus, we can think of L also as a Darboux 
transformation from the operator Laq^b,c,d at the points {xj-^ ,Xj^,..., Xj^ , a+ 
a~^}. From (|4.4p . (|4.5|) . and (|6.2() one can deduce that 

LQ' = Q'Lag,b,c,d with Q' = QQa,b,c,d- (6.24) 

Let us define functions ip'j^i'y) for k = 1, 2, . . . , m + 1 by 

i^kh) = (1 - a/zji^)fj,kr^{aq,b,c,d; Zj^) + i^kS^{aq,b,c,d; zj^) for 1 < A; < m, 

(6.25) 

where //^ and are the constants in (|6.17() . and 

i^'m+iil) = Sj{aq,b,c,d;a). (6.26) 
Using now equations (|4.6a() - (|4.6bl) . we see that 

QaM4^'kil) = ^^^^^Kl^^V'fc(7) for 1 < < m, (6.27) 

which combined with (|6.24|) gives 

ker Q' = Span { (7) , V2 (t) , • • • , C+i (7) } • (6-28) 

Thus, we can think of L as a Darboux transformation from Laq^b,c,d, and 
the kernel of the corresponding intertwining operator is given explicitly by 
(jlT^ and K7^ . 
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Proof of Theorem I6'.M From Remark 16.41 it follows that we can assume m 
is even, i.e. m = 2m for some m E N. Remark 16.31 tells us that without any 
restriction, we can fix the support of Q to be [—rh, +m] (i.e. 7712 = —mi = 
m) in (|6.9j) . and we can choose an appropriate normalization factor 0(7) in 
(|6.9j) . Let us define 0(7) as 

f 1 if ?fi is even; 

Notice that 

I(a(7)) = (-l)'^a(7), (6.30) 
i.e. up to a sign, 0(7) is an /-invariant function. From Lemma l5.2( (|6.17|) . 
and (|6.9j) it follows that Q is a difference operator with coefficients, which 
are rational functions in q'^ , i.e. Q £ TZ{q"',E}. From Lemma [5. 21 and (|6.17|) 
we know that the kernel of Q consists of /-invariant functions, i.e. 

Hi^kil)) = ^kil) for 1 < A; < m. 

From this relation and the definition of the involution /, one can easily check 
that 

I{i^k{l + i)) = ipkil - i) for any i £ Z. (6.31) 
Formulas (|6.9j) . (|6.3flj) . (|6.31j) show that Q is an /-invariant operator, i.e. 
Q G A^. Since Lq = La^b,c,d £ A^, we see immediately from (|6.10p that 
£ G A^, which combined with (|6.1H) gives that V G A^. The proof now 
follows from Theorem 13.31 and Theorem 13.11 □ 

Remark 6.5. Let us denote 

hail) = a'iabq^, acq^ , adq^; q)oo, (6.32) 

and define a new second-order difference operator 2, by conjugating the 
Askey- Wilson operator La^b,c,d in (!2.1j) with /ia(7), i-e. 

Sl = ha{7r' La,b,c,dha{l)- (6.33) 

If we write £, in the form 

£ = + *B^Id + (t^E-^, (6.34) 

then a straightforward computation shows that the coefficients of £ are given 
by the following formulas 

_M7+1) 1 - abcdq'-' 

^~ hai-f) ^ ~ (1 - abcdq^"f-^){l - abcdq^^) ^^'^^^^ 



^ = Ml^lc^ = (1 _ qi){i _ abq^-^){l - acq"<-^){l - adq^'^) 

hah) 

{l-bcq^~^){l-bdq-y-^){l-cdq^-^) 

(1 - abcdq^^-^){l - abcdq'^y-^) ^ ^ ^ 

_ (1 + abcdq^^-^){sq + s'abcd) - g^"^(l + q)abcd{s + s'q) 

}5^-B^-q {I - abcdq^"i-^){l - abcdq^^) 

(6.35c) 
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where s = a+b+c+d and s' = a^^ +b^^ +c^^ +d~^ . Notice that the operator 
£ is symmetric in o, b, c, d (in fact, this was the operator originally introduced 
by Askey and Wilson, see 2, p. 5]). Therefore, we can produce bispectral 
operators as in Theorem 16.21 if any two of the parameters (a, b, c, d) are 
related in the same manner as a and d in (|6.15|) . 

Theorem 16 . 21 and E,emark l6.5l tell us that if two of the parameters a, b, c, d 
are related by a formula similar to (|6.15|) . we can construct bispectral op- 
erators by adding bound states to the left or to the right of the continu- 
ous spectrum, depending on the sign of these two parameters (they must 
have the same signs!). Assume now that we have two positive parameters, 
satisfying 1)6. 15() with -|-, and the other two parameters are negative and 
satisfy (|6.15|) with — . Then it is possible to construct bispectral operators 
by adding bound states both to the left and to the right of the continuous 
spectrum. The corollary below explains what modifications are necessary 
in the proof of Theorem 16.21 in this case. The operators considered below 
contain the so called continuous g-Jacobi operators which can be obtained 
from the Askey- Wilson's operator by taking 

a = q"^^, b=—q^^^, c=—q^, d = . 

This substitution was used by Rahman in For a different parametriza- 
tion of the continuous g-Jacobi polynomials, as well as relations between the 
different parametrizations see (23 p. 83]. 

Corollary 6.6. Assume that 

a = qh/2+a^ 6=_gW2+/3^ c = -qh/^ ^ d = qh/2^ (6.36) 

where a,(3,li,l2 are positive integers. Let us define 

J = |g«i/2+fc + g-'i/2-fc, < /fc < a - l}, 

K = I - (g'2/2+fe + q-l^l^-ky < < /? - l}, 

and let X = {xj^ ,Xj^,..., Xj^ } be a subset of JUK. The second-order differ- 
ence operator L which is obtained by iteration of the Darboux transformation 
1)6. If) , starting with the associated Askey-Wilson operator Lq = La^b,c,d (|2.1|) . 
and adding bound states at the points of X possesses a two-dimensional space 
of eigenf unctions which are also eigenf unctions of a q-difference operator in 
the spectral variable z. 

Proof. As in the proof of Theorem 16.21 we can assume that m is even, i.e. 
m = 2m and the support of Q is [—rh,rh]. From Lemma 15.21 we know 
that the eigenfunctions of L, corresponding to eigenvalues Xj^. £ J are I- 
invariant. The main point is to see what happens with the eigenfunctions 
■0fc(7), corresponding to Xj^, G K. Let us denote by £^ the Askey-Wilson 
operator with parameters a, b, c, d replaced by b, c, d, a, respectively, i.e. 



(6.37) 
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Then using the notations in E,emark 16.51 we will have 

£ = /ife(7)-^£'/ife(7), (6.38) 

where /ife(7) is analogous to ha{'j) defined in (|6.32|) : 

h{j) = P{baq\ hcq\ hdq^; q)^. (6.39) 

From (IOhI) and (lOsl) it follows that 

J _ hail) b hbil) I 
^a,b,c,d - -r-r^ ^ , / V (b.4Uj 

hh) hail) 

This equation simply means that the correspondence between the eigenfunc- 
tions of £^ and Lafi^c,d is given by multiplication by the function /iq (7) 7/1^(7), 
i.e. if 

then 

ia,fe,c,d^a(7) = A^'a(7) if and only if £^^6(7) = X^bil)- 

If we take an eigenfunction of with eigenvalue xj^, G K, then this function 
must be /-invariant by Lemma [5. 21 interchanging the roles of the parameters 
ia,d) and (6, c). Thus, it remains to see how /i^ (7) //ife (7) is transformed 
under the involution /. Using (|6.3()j) . (|6.32j) and (|6.39j) we get 

Mt) ^ (ayiacqTMqr^ (-1)^/^(7) (6 41) 

hil) \b) ibcq^hdq^-q)^ ^ ^^^>' ^^"^'^ 

where 



(g7+/3+«2^ _q7+(«l+«2)/2+/3. 



00 



Using the definitions at the beginning of Section 2 we can write /i(7) in the 
form 

q[{h-l2)/2+a-fS]-y 

^^^^ ^ (g^+^+'^ ; q)a-,+l.-l, (-,7+fc+^2)/2+/5. ,)^_^ ' (6-43) 

where 

^^^'^)"=(^^'°^"<°- 
Notice that if li — I2 is even, then hii) is a rational function in q^ , i.e. 
/i(7) S TZ{q'^,E}. However, for li — I2 odd, /i(7) has an extra factor q^^"^, 
i.e. /i(7) = g"^/^ X (rational function in q^). 

In the case considered here (i.e. assuming (|6.36p ). the involution / acts 
on rational functions in q^ by the following formula 

/(q7~) — q-^-a-P-h-h+'L _ 

We can naturally extend the involution I to an involution acting on TZ{q'^/'^^ 
by defining 
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Using (|(i.4:-{|) and (|(i.44|) one can check that 

= (6.45) 

i.e. up to a sign the function h{'y) is /-invariant. Finally, let us denote by m' 
and m" the number of elements in the sets X f] J and X D K, respectively. 
The computations above can be summarized as follows: 

• m' of the functions ipk, generating the kernel of Q, are eigenfunctions 
of La,b,c4 corresponding to eigenvalues Xj^ G J and these functions 
are /-invariant rational functions in q^] 

• m" of the functions 'i/'fc, generating the kernel of Q, are eigenfunctions 
of La,b,c4 corresponding to eigenvalues Xjj_ G K and these functions 
are of the form 

Ml) = {-^Vikii), 

where 

mm = {-ir^''^''^''Ui), (6.46) 

and 

, , _ J rational function in if li — I2 is even; , , 

= I qil2 X (rational function in q'l) if h - I2 is odd. ^ ' 

Without any restriction, we can assume that Xj^ £ J for k = 1,2, ... ,m' 
and Xj^ £ K for s = m! + l,...,m". Let us denote by {aij;hij;ci) the 
(m -|- 1) X (m + 1) matrix with rows numbered from rh to —in (i.e. i = 
fh,fh — 1, . . . , —rh) and columns from 1 to m -|- 1, with entries aij for the 
first m' columns, entries hij for the next m" columns, and entries Cj in the 
last column. Then, from (|6.9|) we see that the operator Q has the form 

Q = a(7)(-l)"^"^ det(^,(7 + i); (-1)^0(7 + »); E'), (6-48) 
where 0(7) is an appropriate function. From ()6.46|) it follows that 

which combined with (|6.4<S|) shows that 

/(Q) = (-l)^^/(a(7)(-l)™"^)det(V^,(7 + i);(-l)'e,(7 + i);i?'), (6.49) 
where M = m + {a + l3 + li + l2)m". From (|07|) and we see that 

_ , ^Nm"7 / operator in T^jg''', i?} if m"(/i — /2) is even; 

-2 = a[7){-l) X I ^^/2 X (operator in n{q-^,E}) if m"(/i - I2) is odd. 

(6.50) 

Let us now define 0(7) as follows 

" (_l)m"7(g^l7/2 + g2g"'i('^+"+'^+'i+'^"^)/^) ' ^^'^^'^ 

where 

2 if m"{li — I2) is even; 



ei = 

and 



1 if m"{li — I2) is odd; 

^2 = (-1) 



1 if M is even; 
1 if M is odd. 
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From (|OH|) . (jnSSl), (|H3n|) it follows that with this definition of a{j), Q 
becomes an /-invariant operator from TZ{q'^,E}, and the proof continues as 
in Theorem 16.21 □ 



7. A (^-ANALOGUE OF THE KrALL-JACOBI POLYNOMIALS, WITH A BOUND 
STATE OFF THE CONTINUOUS SPECTRUM 

In this section we illustrate Theorem 15. 21 on the special case 

a = q2 , d = , 6 and c arbitrary, (7.1) 
adding a bound state off the continuous spectrum of La^b,c,d at the point 

xo = q^+q~^. (7.2) 
This corresponds to the choice 

1 = 1, a = l, k = 0. (7.3) 

7.1. Illustrating the general theory. According to Remark 16. 41 this case 
must be handled by first expressing La^b,c,d as a Darboux transform (without 
free parameter) from Laq^b,c,d at a + a~^, and then performing a Darboux 
transform (with free parameter) of La^b,c,d at ^o- Thus we perform the chain 
of elementary Darboux transformations 

Laq,b,c,d = {a + a ^) Id + Pa,b,c,dQa,b,c,d ^ 

La,b,c,d = (a + a~^) Id + Qa,b,c,dPa,b,c4 = xo Id + PQ 

rxL = XQld + QP. (7.4) 

We now explain how Theorems 3.1 and 3.3 (which form the core of the 
general proof of Theorem 16. 2|) apply to this example. We can factorize 

^ = {Laq,b,c4 - 3;oId) {Laq,b,c,d - (o + a~^)Id) , (7.5) 

as 

C = VQ with r = Pa,b,c,dP and Q = QQa,b,c,d- (7.6) 

The explicit formulas for Pa,b,c,d and Qa,b,c,d are given in (|4.2() and 1)4. 3(1 
respectively. The kernel ip^ of the operator La^b,c,d — xo Id is readily obtained 

by substitution of (f7?T|) and (f7^ in ((^ and ((^ (using (jSHJ), (f^ 
and 1)5. 7|l '). which leads to 

"^'^ {I- bq^+'){l - cq-^^l) [^^{l- lAil - bcq^))' ^^'^^ 
with /u a free parameter. This gives 

Q = V.,(V^,i?-^,+i), p = _9j.^E-^ + ^:p^. (7.8) 

In this last formula is as in (|4.H) and corresponds to the factor 0(7) 
defined in (|6.29|) (with fh = 1 and a replaced by aq), and A^.a,b,c,d and 
Cry-a,b,c,d are defined as in (|2.2aj) and ()2.2b|) : it is understood that the values 
(|7.1j) for a and d are to be substituted everywhere. 
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The following identities are equivalent to (|4.5|) : 

-!fj— A — A nnrl Jfj—O —C 



^'•f+1; a,b,c,d ■^'y; aq,b,c,d anu ^'y; a,b,c,d ^'y; ag,b,c,d- 

Using these identities, by a straightforward computation, we obtain that 

Q = QQa,b,c,d = A^;aq,b,c,d '^7 E 

— {Ar,a,b,c,d '07+1 + Cf+l; a,b,c,d i^-y) W + C^-aq,b,c,d '07+1 E'-^ , (7.9) 

and 

V = Pa,b,c4P = 

( ^7;o,b,c,d _^ a,b,c4 j _|_ ^7; aq,b,c,d 

Ip^ V7+I / ^^7 

Since C in (|7.5() is a constant coefficient polynomial in Laq^b,c,d, the relevant 
involution / is obtained by replacing a with aq in H3.15() . that is 

I{qr) = -J— = and I{E) = E'K 

^ ' abcdq-y bcq-y+^ ^ ' 



We know from (|3.18|) that l[Ay^aq,b,c,d) = C^;aq,b,c,d- Using (|2.2aj) . 
and (|7.7|) . it is straightforward to check that 

I{A^;a,b,c,d) = C^+l;a,b,c,d and = V'7+l ■ (7-11) 

This makes clear that the operators Q and V as expressed in (|7.9j) and (|7.10|) 
are /-invariant. To alleviate the notations, we shall write L and A to denote 
the following difference operators: 

L = Lag^b,c,d, ^ = — q'^){^ ~ abcdq'^) Id, with a,b,c,d as in ()7.1() . 

According to Theorem 13.31 we must have that 

r = ur-^ and Q = e-^v, 

with U,V eB = {L, A) and T, 6 G /C = (A). 

Using the symbolic software MACSYMA, we have obtained the following 
explicit formulas for T,Q,U and V. It follows from (|7.11|) that the function 

"07+1 is /-invariant, hence it can be written as a rational function of A. 
Explicitly: 

^P^ipy+i Id = r(A)e-^(A), with 

r(A) = [{q - l){bc - l)/i + 1] [{q' - l){bcq - l)/i + g]ld 
+fi[q + l + fi{q - l){2bcq - q + bc-2)]A + IS? , 
e(A) = q-^ [A + (g2 - l)(6c - l)Id] [A + (g - l){bcq - l)Id] 

X [A + (6gl/2 _ l)(cg3/2 _ ^ (^^3/2 _ l)(cgl/2 _ 1)1^] _ 
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Then, one finds that 

U = {riK + r2h? + n A^) L + L (si A + sa + S3 A^) 

+ t L + uqM + ui A + U2 + U3 A^, 
y = (wi A + t;2 A^ + ^3 A^) L + L {wi A + W2 A^ + wz A^) 

+ x A L A + 2/0 Id + yi A + A^ + yg A^ 

for some appropriate choice of the constants rt, Si,t,Ui,Vi,Wi, x,yi. The 
exphcit expressions for these constants are rather comphcated and there is 
no point in displaying them. Theorem 13. II fsee especially (|3.1U|) and (|3.11|l ) 
provides the following explicit bispectral operator 

B = b{V)b{U)r\ (7.12) 

with 

/ = b(£) = {z + z-^ - {q^ + q-^)) {z + z-^ - {q'l + q'^)), 
b{U) = iz + z-'^){nB + r2B^ + rsB^) + (siB + ss^^ + S3B^)iz + z"!) 

+ t{z + + Mo + UiB + U2B'^ + u-sB^, 

b{V) = {z + z~^){viB + V2B'^ + v^B^) + {wiB + W2B'^ + w-iB^){z + z"^) 

+ xB{z + z-'^)B + yo + yiB + ya^^ + y^B^ , 

where B = -6^5/2 ;, c,gi/2(-S5 D^) and Ba^b,c,d{z, Dz) is defined as in ()2.6p . Since 
-B is a second-order (^-difference operator, B in (|7.12p is a (/-difference oper- 
ator of order 12. 

7.2. A bispectral operator of order 4. It is possible for the example 
discussed above, with a, b, c, d as in 1)7. If) , to produce a bispectral operator of 
order 4, by performing a bispectral Darboux transformation of the operator 

>c = (V^'>,c,,v2-(g'/' + g-'/')id)', 

instead of the choice we made in 1)7. 5|) . The idea is to first obtain the operator 
Lg3/2ftcgi/2 as a (contiguous) Darboux transform of ^ f. gi/2 at q^^'^ + 
using ()4.4|) and 1)4. 5|) . and then to perform a (general) upper- lower 
Darboux transform (with a free parameter) of ^^3/2 at q^^"^ + 

In the case q = 1, such a strategy has already been explained with all details 
in Section 4 of ^H]- Since we don't know yet how to produce a bispectral 
operator of the lowest possible order for the general situation described in 
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Theorem \H.'2\ we shall only display the final result. Let us define 

/i q^^z -b){z- c){z - bq){z - cq){z - q^l'^){(i^l'^ z - 1) 

~ {Z + ^l/2)(^ + g3/2)(^2 _ l)(^2 _ ^2) - 

D-i{z) = q{q + l){z - b){z - c){z - q^'^) 
[z + q^/^){q''^z + 1) - \{q + bc){z^ + q) - (b + c){q^ + l)z 



{z + q^/^){z + g3/2)(gl/22 + i)(^2 _ 1) 
=D_i(l/z), D2{z)=D^2{l/z), 
Do{z) = -D_2{z) - D_i{z) - Di{z) - D2iz), (7.13) 

and 

(g^+l - l)(6cg^+l - 1) \q^'{q + I) + fi {bcq< - - 1) 

A(7) = ^ 
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Then, the functions 

\E'(7,z) = (5^'(7,2:), with 

^(7,z) = i?^(g3/2,6,c,gi/2.z) or 5^(g3/2, ft, c, g^/^. z), 

and with Q as in 1)7. 8() . R-y{a,b,c,d; z) and S-y{a,b,c,d; z) as in 1)2. 8a() and 
(fT8b|l . satisfy the pair of equations 

L^(7,2;) = (z + l/z)l'(7,z), Las in ((EH), 

2 

J] A(^)^'(7,9'^) = A(7)^'(7,^)- 

i=-2 

Let us consider the special case of (|7.1() given by 

a = q'/^, d = q^'\ b = -q^l\ c = -q^''^'^ (5 > -I, (7.14) 

and let us also specialize 7 to be a positive integer. Then, the functions 

R^{q^'\-q^'\-q^^''\q"^-zl n = 0,1,2,..., 

become (up to a normalizing factor) the continuous g-Jacobi polynomi- 
als with parameters (a = 1, /3) as defined by Rahman (see [25!) formula 
(3.10.14)). When g — > 1, these polynomials reduce to the standard Jacobi 

polynomials Pn'^\x), x = {z + z^^)/2 (see [25^, formula (5.10.2)), with 
weight function (1 — x)(l + x)^ on [—1, 1]. The new functions 

Pg{z) = 1, 

Pn{z) = {QR)n-l = ^n-li^n-lRn{q^^\-q^^^-q^+^^\q'^^;z) 

-0„_iV^„i?„_i ((7=^/2, -gV2, _g/3+i/2^ qi/2. n = 1, 2, 3, . . . , 

obtained from the Darboux process, with Q as in (|7.8)) . are (Laurent) poly- 
nomials of degree n in the variable z+z~^. These polynomials are eigenfunc- 
tions of the fourth-order g-difference operator with coefficients as in (|7.13|) , 
with b and c as in (17.141). 
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Putting the free parameter /i of the Darboux transform in (|7.7|1 to be 



^ 2/3+1(^-1)2' 
one can check that, taking the hmit q ^ 1, the polynomials 

z -\- 

Kn{x) = Hm (g - l)^P„(z), x = , 

are the so-called Krall-Jacobi polynomials, which are orthogonal on the in- 
terval [—1,1] for the weight function 

{l + xf + p5{x-l), (7.15) 

where 6{x) denotes Dirac's delta function. This is one of the family of or- 
thogonal polynomials satisfying a fourth-order differential equation that was 
discovered by H.L. Krall in j^H]- Thus, the example discussed in this section 
can be viewed as a g-deformation of the Krall-Jacobi polynomials, with a 
bound state of arbitrary weight p off the continuous spectrum, precisely at 
the point x = {z + z~^)/2 = {q^/"^ + q~'^/'^)/2. When g ^ 1, this bound 
state tends to the boundary x = 1 of the continuous spectrum of the Jacobi 
polynomials, in agreement with (|7.15)1 . 

8. AsKEY- Wilson type solitons 

In this section we construct a large family of second-order difference op- 
erators which provide solutions to the Askey- Wilson bispectral problem as 
defined in (|1.1|) and (|1.2j) . within the context of the theory of commutative 
rank 1 rings of difference operators. We start from the simple observation 
that when a = — 6 = 1 and c = —d = y^, the Askey- Wilson operator (|2.1|) 
reduces to 

La,b,c,d = E + . (8.1) 
More precisely, we shall prove the next theorem: 

Theorem 8.1. Let ki,l < i < g, denote g positive integers. Let us par- 
tition these integers into two arbitrary sets J and K of a and (3 elements 
respectively, a, /? > 0, a -|- /? = We build accordingly the a points q^^/'^ + 
E J, and the /3 other points —{q^^/'^ + q~^^/'^),i G K, and we assume 
that all these points are distinct. 

i) The tridiagonal operator L which is obtained by iteration of the Darboux 
process starting at L = E + E~^ , and adding bound states at these g = a + j3 
distinct points, is part of a commutative rank 1 ring of difference operators, 
with spectral curve as in (|1.6|) . 

a) The (reduced) wave function ^(n,z) of this rank 1 commutative ring 
of difference operators, is also a common eigenfunction of a rank 1 commu- 
tative ring of q- difference operators in the spectral variable z. 

The proof of this result is closely connected with the idea of " duality" that 
we exploited to establish Theorem l2.lL In the context of rank 1 commutative 
rings of differential operators, "duality" was first systematically exploited 
by G. Wilson Before proving Theorem 18. H we need to give a brief 
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summary of the theory of rank 1 commutative rings of difference operators. 
We Umit ourselves to the case when the spectrum of the ring is a rational 
irreducible curve, which is the only part of the theory that we shall use. The 
presentation follows closely our previous work jlBj . 

8.1. Commutative rings of difference operators and rational curves. 

We denote by "P = C[z] the space of polynomials in the variable z, and we 
denote by 7^ = C{z) the space of rational functions of z. 

Definition 8.2. A flag V : • • • C Vn+i C Vn C Vn-i C • • • of subspaces 
Vn C TZ is called a rational flag if and only if there are polynomials p{z) and 
r(z) (independent of n) such that 

z'^p{z)V CVnC z'^r-\z)V, 

and the codimension ofVn in z"'r~^{z)T' is equal to the degree ofr{z). 

There is a one-to-one correspondence between rational flags and affine 
irreducible rational curves which complete by adding two non-singular points 
at infinity. The curve is recovered from the flag as the spectrum of the ring 

Ay = {f{z) £ TZ with poles only at z = and z = oo, 

such that 3k £ Z for which f{z)Vn C Vn}. (8.2) 

When the curve is rationally parametrized by z, the two non-singular dis- 
tinguished points at infinity Poo and Qoo correspond respectively to z = oo 
and z = 0. 

By definition, the rank of a commutative ring of difference operators is 
the greatest common divisor of the orders of all the operators in the ring. 
There is a dictionary between the affine rings of irreducible algebraic curves 
which complete by adding two non-singular points at infinity and the rank 
1 commutative rings of difference operators. We refer the reader to |31j for 
a full account of the theory, as well as for adequate references. We limit 
ourselves to the case of irreducible rational curves as described above. The 
general case proceeds along similar lines. To each Vn there corresponds a tau 
function ry„(t) = T{n,t) in the sense of Segal and Wilson 36 , depending 
on infinitely many variables t = {ti,t2,t3, . . .). The wave and adjoint wave 
functions, denoted respectively by w{n,t,z) and w*{n,t,z), are defined by 

win, t, z) = z"exp(^, z) \~ (8.3) 

T[n, t) 

w* (n, t,z)= z-"exp-i (t, z) ^("'^ + [^ ^]) ^ ^j^h (8.4) 

T{n,t) 

exp(t, z) = exp(tiz + t2z'^ + t^z^ H ) and [z] = [z, z^ /2, z^/3, . . . ). 

The plane Vn can be recovered from the wave function as 

Vn = Span{t(;(n, 0, z),w{n + 1, 0, z),w{n + 2, 0, z), . . . }. 
The dual flag of subspaces V* : • • • D V*j^^ D V* D V*_-^ ^ • • • , with 
V* = Span{w*(n, 0, z),'w*{n — 1,0, z),w* {n — 2,0, z), .. . }, 
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defines the same tau function, up to a change of sign, Tv*{t) = Ty„(— t). 

Theorem 8.3. (see ^H], and references therein). The wave function 
w{n,t,z) is the common eigenfunction of a commutative rank 1 ring of dif- 
ference operators which is isomorphic to A^, as defined in ()8.2|) . More 
precisely, for any f € A^, there is a difference operator Lj such that 

Lfw{n,t, z) = f{z)w{n,t,z). 

If f has a pole of order i at P^o and a pole of order j atQ^o, the operator Lj, 
thought of as a finite hand matrix, has i diagonals above the main diagonal 
and j diagonals below it. 

We illustrate the concepts above on the example of the simplest rational 
singular curves, namely those which are obtained by identifying 2g distinct 
points on the Riemann sphere P^(C) in pairs {Ai,/Xj},i = \,...,g. We 
assume that all these points are distinct from the two distinguished points 
Qoo = and Pqo = oo. We pick g nonzero arbitrary complex numbers 
(5i G C*, 1 < i < 5 (i.e. a divisor of degree g on the curve), and we define 

Vn = :p=fo — z --[meromorphic functions / on P"'^(C) such that 

(/)-riQoo >OonPi(C)\{oo} and /(A,) = 5J(/Xi), i = l,...,5}. (8.5) 

By the notation (/) — nQoo > 0, we mean that / has a zero at least of order 
n at if n > and a pole of order at most —n at Q^o if < 0. The 
next lemma shows that this definition fits within Definition 18.21 In order 
to establish it, it is useful to introduce the following definitions. A function 
/ G 7^ which admits an expansion around z = cxd of the form 

/ = c„z" + c„_iz"-i + • • • , c„ / 0, (8.6) 

will be called an element of order n. It will also be useful to consider the 
non-degenerate bilinear form on IZ 

B{f, g) = ms,=^f{z)g{z)dz, f,g£n. (8.7) 
Lemma 8.4. The plane Vn defined in 1)8. 5(1 satisfies 



a 

z 

1=1 



" ri(^ - c K c \{{z - \i)-'v, 



and the codimension ofVn in z"'Y\f=i{z — Aj) is equal to g. 

Proof. The second inclusion as well as the assertion about the codimension 
both follow immediately from the definition of Vn in (|8.5)) . In order to 
establish the first inclusion, we introduce the dual flag of subspaces 

V* = — ^ ^|meromorphic functions / on P"'^(C) such that 

(/) + nQoo > on \ {oo} and fipa) = 5*/(A,), i = 1, . . . (8.8) 
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with 

It follows immediately from (|53|) . and that 

{iesz=\i + 'i:(iSz=fjLi]w{z)w* {z)dz = 0, V G a.nd'i w*{z) £ V^. 

For G l/„ and G V^^ with n > m, w{z)w*{z) has no pole at 2; = 0. 

Thus, by the residue theorem, we deduce that 

9 

r:eS;,=oow{z)'w*{z)dz = - y^{res^=A, + TeSz=f^Jw{z)'w* {z)dz = 0, 

i=l 

V w{z) £ Vn and V ■w*{z) £ with n > m, 

showing that V* C V^, with the orthogonal of Vn with respect to the 
bilinear form B defined in H8.7() . Since has a basis of elements of orders 
n, n+1, n+2, . . . (see (|8.6() for the definition), and V* has a basis of elements 
of orders — n,— n+2,...,we deduce that V* = V^. By the definition 
of V*, it imphes that 

9 9 

= V: C z-^ J{{z - f,,)-'V ^VnDz^ fliz - f,,)V, 

i=l i=l 

as desired. This finishes the proof of the lemma. □ 

Lemma 8.5. The tau function T{n,t) associated with the rational flag V 
defined by (|8.5|) . is given by 

9 00 g 

rin, = n /^7"e^P ( - E E ^0 

j=l i=l j=l 

dei(\'l^'-^exp{t,Xj) -6jn''+'-^exp{t,^ij)) . (8.10) 

Proof. It follows easily from the definition of Vn in (|8.5p that the wave 
function associated with this plane is 



exp(i,2;) 

w[n, t, z; — 



nj=i(^-A,) 

det(A^"+^-^exp(t, A,) - 6jfi]+'~^exp{t, z^+'-^) 
det(Aj+*-^exp(t, A,) - 6,^i]+'-'expit, ^i,)) ^^^^^^^ 



i.ll) 



The notation in the numerator means the determinant of the (^ + 1) x (g + l) 
matrix with {i,j) entries as indicated there for 1 < i < g + 1,1 < j < g, 
and entries z"+*~^,l < i < 5 + 1, in the last column. Elementary row 
manipulations with this determinant lead to 

_ eMt,z) ^n-,, nU^-^Mn,t-[z-^]) 
"'^"'''"^-n?=i(--A,) .Mn,t) 
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in agreement with (|8.l-{)) . which defines the tau function. This finishes the 
proof. □ 

Remark 8.6. The motivation for the (irrelevant) factor 11^=1 ^7" the 
definition 1)8. 1U() of r(n,t) wih appear in the next subsection, see (|8.36l) . 

Having in view Theorem 18.11 we are particularly interested in rational 
curves with double points, for which the associated commutative rings of 
difference operators contain a tridiagonal operator. 

Proposition 8.7. Let us assume that fii = A~^, and that the 2g points 
'^ii/^ijl ^ i ^ g, o-fG still distinct of each other, i.e. Aj ^ ±1 and Aj 7^ 
A^"'^,V i ^ j. Then, the rank 1 commutative ring of difference operators 

Ay, with Vn as in ()8.5() . contains a tridiagonal operator L with one diagonal 
above and one diagonal below the main diagonal. The spectral curve of the 
ring has for equation 

9 

Spec(^v): = {u' -A)J[{u-{\i + \T^)f, (8.12) 

i=l 

and the operator L can be obtained by iterating the Darboux process starting 
with L = E + as in (|8.1|1 . and adding bound states at the g points 
Xi + X-\ 

Proof. Since Hi = A~^, it is clear that the function 

ii = z + z-i, (8.13) 

belongs to the ring Ay as defined in ()8.2() . with uVn C Vn-i- This function 
has a simple pole both at z = and at 2: = 00 hence, by Theorem l8.31 there 
exists a tridiagonal operator L„ with one diagonal above and one diagonal 
below the main diagonal, satisfying Luw{n, t, z) = uw{n, t, z). 

Since u has a simple pole at z = 0, all other generators of Ay can be 
taken to be polynomials. The functions 

9 

fk = z-'^X[{z-\i){z-\r^), keZ, (8.14) 

i=l 

which vanish at A, and /ij = A~^, obviously belong to Ay. For I < k < g — 1, 
by subtracting an appropriate polynomial of degree k in u, one obtains 
polynomials (in z) q2g-i{z), q2g-2{z), ■ ■ ■ , qg+i{z), of degrees 2g - l,2g - 
2, . . . , g+1, that belong to Ay. Combining them with the functions fk, k <0, 
defined in (|8.14j) . we conclude that Ay contains polynomials qk{z) of degree 
k, \/k > g -\- 1. On the other hand, any polynomial q{z) £ Ay of degree 
k < g, must be identically constant. Indeed, such a polynomial must satisfy 
q{Xi) = g(A^^),V 1 < i < g, which amounts to a linear homogeneous system 
Ac = 0, for the unknown coefficients ci, . . . , Cg of q{z) = Co + ciz + - ■ ■ + CgZ^. 
One shows easily that 

9 9 

detyl = nAr^n(Af-l) \{ (A, - A,) \{ (1 - A,A,). 
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Since the 2g points Aj,;Uj = Aj" , 1 < « < 5, are assumed to be distinct, this 
determinant is nonzero, implying that ci = C2 = ••• = % = 0. 

The upshot of the discussion above is that the algebra A\; is generated 
by u in (|8.13j) and polynomials of degree k, k > g + 1, ot equivalently by u 
and the functions fk,k G Z, introduced in (|8.14j) . We now show that the 
functions u and 

i=l 

are enough to generate Ay, which will establish ()8.12|) . Clearly 

9 
i=l 

Since ufk = fk-i+fk+i, we deduce inductively that fg-i, fg+i, fg^2, fg+2, 
belong to the algebra generated by u and v. 

It is well known that the spectral meaning of the curve H8.12() is that the 
operator = L, tt as in H8.13() . can be obtained by iteration of the Darboux 
process, starting from L = E + and adding bound states at the points 
Xi + A-^ This can be easily deduced from the explicit form for the wave 
function given in (|8.11|) . with /ij = A^"^. The proof is complete. □ 

8.2. The proof of Theorem 18. IL The proof of Theorem 18.11 is a combi- 
nation of Proposition 18.71 above and Proposition 18.81 below. It is convenient 
to introduce the so-called reduced wave and adjoint wave functions, which 
are obtained by omitting the factors exp(t,z) and exp~^{t,z) in 1)8. 3() and 
1)8. 4() respectively: 

^{n, z) = ex.p~^{t, z)w{n,t, z), ^*{n, z) = ex]i{t, z)w*{n,t, z). (8.15) 

We also omit to write the explicit dependence on t, which is irrelevant for 
what follows; we just think of ti, t2) ^3; ... as some free parameters. 

Proposition 8.8. Consider a rational curve with double points, obtained by 
identifying 2g distinct points of¥^{C) in pairs {Aj,//j},i = 1,... ,g. Assume 
that 

Xi = q^-^ii, with /ci G N \ {0}. (8.16) 

Then the reduced wave function ^{n,z), besides being a common eigenfunc- 
tion of a rank 1 commutative ring Ay of difference operators (in n) 

Lf^{n, z) = f{z)^{n, z), V/ E Ay, 

is also a common eigenfunction of a rank 1 commutative ring of q- difference 
operators A<^, in the variable z 

Bg^{n,z)=g{n)^{n,z), G Ay. (8.17) 

Assuming for a moment this result we can prove Theorem 18. II 

Proof of Theorem \8.1\ From (|8.16() we immediately obtain that the only in- 
tersection between Proposition 18.71 and Proposition 18.81 is when 

Xi = q'^'m and m = X^^ ^ Xi = ig'''/^ /cj G N \ {0}, 
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in which case (|8.12j) reduces to □ 

As already announced, the proof of Proposition EHl exploits an interesting 
duality of the (reduced) wave function. 

Proof of Proposition \8.8[ Using (|8.16() , we have that 

fi^lx^^X'r^expit, Xj) - (5j^}-^exp(t, /ij)}, with x = g". (8.18) 

From (|8.11|) it follows that the (reduced) wave function is a function of 
X = and z, i.e. 



*(n,x) = ^'(x,2;) = 

^iog_^i2^ clet(x'^J A*~"'^exp(t, Xj) - 5jf/j'^exp{t, Hj); z^~^) 
nf=i(^ - Ai) det(x'=^ A;.-iexp(i, A,) - <5,/i;.-iexp(t,//,))^<. 
Any (reduced) wave function admits an expansion 



.19) 



log X log z / Wlix) Woix) \ 

^(x,z)=e 9 1 + — ^ + ^V^ + --- asz^oo. (8.20) 

\ Z Z"^ / 

Since the ki are positive integers, it follows immediately from (|8.19|) that 

log X log z det(A!,~^: z*"^) 

lime ^^^^^^(x,z) = — ^ = 1, (8.21) 

nf=i(^-A,)det(A;.-i)i<,,,<3 

showing that, when the special conditions ()8.16p are satisfied, we also have 

N log X log z f wAz) W2(z) \ 

^{x,z) = e i°g 9 1 + ^ ^ + , H asx^oo. 

This suggests that the function ^'(x, z) obtained by exchanging the variables 
X and z in ^'(x, z) 

^{x,z) = ^{z,x), (8.22) 

should be the (reduced) wave function associated with another rational flag 
V : • • • C Vn+i C Vn C Vn-i C • • • , in the sense of Definition 18.21 The cor- 
rectness of this assertion will be established in Lemma l8 . 91 b elow . Assuming 
the result, we deduce from Theorem 18^ that 



Sg§(x,z)= Di{x)i!{q'x,z)=g{z)i!{x,z), g{z) e A^, 

finitely many j's G Z 

which, remembering the definitions of ^{x^z) and ^{x,z) in ()8.19() and 



(|8.22|) . amounts to 

^ D,{z)^{n,q'z)=g{q^)^>{n,z), 

finitely many i's G Z 

which establishes (|8.17|1 . This finishes the proof of Proposition 18.81 □ 
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As mentioned in the proof of Proposition 18.81 given above, to complete 
the argument, we stih need to estabhsh that ^{x,z) defined as in (|8.22j) is 
indeed the (reduced) wave function of some rational flag. We need some 
preliminaries. 

We introduce the following multiplicative groups of formal pseudodiffer- 
ence operators 

oo oo 

W = {l + ^ Wijx-W-'Y W* = [l+Y, WijX~Wi'}, (8.23) 

i,j=l i,j=l 

with Wij G C, Dx and D^^ the forward and backward g-shift operators, 
Dxh{x) = h{qx) and Dx^h{x) = h{q^^x). 
We denote by a the adjoint isomorphism 

a:W ^W* -.W r^a{W) = {W~^)*, (8.24) 

where * denotes the adjoint operator, (D^.)* = -D"*. We define an anti- 
isomorphism b : W — > W by 

b(x) = Dx, b{Dx) = X, i.e. 

oo oo 

b(l + ^ WijX'Wx'^ = l+Y^ WijX-'D-^, (8.25) 

i,j=l i,j=l 

as well as an anti-isomorphism b* : W* W* by 
b*ix) = D-\ b*iDx) = x-\ i.e. 

oo oo 

b*(l+Y^ WijX-Wi^ =1+Y^ WijX-'Di. (8.26) 

i,j=l i,j=l 

With these definitions, it is straightforward to check that 

b*a = ab. (8.27) 

We can write (I8.2U|) as 

log X log z 

^(x,z) = W{x) e i°s 1 , 

W{x) = 1 + wi{x)Dx^ + W2{x)Dx'^ + ■■■ . (8.28) 

The functions Wi{x) are rational functions of x which, remembering (|8.21l) . 
satisfy linix^oo 'Wi{x) = 0. Thus 

oo 

W{x) = 1 + X] 'WijX'^Dx' as X ^ OO, (8.29) 

which shows that W{x) S W as defined in ()8.23() . 

It is a general fact from the theory of the discrete KP-hierarchy (see [TFj , 
[201) the adjoint (reduced) wave function as defined by 1)8. 4() and 1)8. 15() 
can be expressed in terms of W{x) as follows 

log X log z 

'i/*{x,z) = a{W){q^^x) e isf^— , (8.30) 
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with a{W) defined as in (|8.24|) . In the same papers, one can find a proof 
that (|8.28|) and (|8.30|) imply the so-called bilinear identities 

B{^{q'x,z),'^*{x,z)) =0, Vi>0, (8.31) 

where B is the residue pairing introduced in H8.7|) . 

With these preliminaries, we can establish the last needed lemma. 

Lemma 8.9. With the assumptions of Provosition TKR. the new functions 

^{x,z) = -^{z,x) and ^* {x, z) = {qz,q''^x)q~^xz~^ , (8.32) 

are respectively the (reduced) wave and adjoint wave functions associated 
with a rational flag V. 

Proof. In terms of the anti-isomorphism b introduced in 1)8. 25|) . ^{x,z) as 
defined in ()8.32|) can be written as follows 

_ log X log z 

^{x,z) = b{W){x) e '"g-J , (8.33) 

with W{x) as in Using (I071) and the definitions and 

of b* and ^*{x,z) respectively, one checks easily that 

i!*{x,z) = b*{a{W)){q-^x) e~'°'i°g'°/ ' , 
from which, using ()8.27p , we deduce that 

i!*{x,z) = a{b{W)){q'^x) e~'°'i°g'r (8.34) 

Let us define 

Vn = Span{§(n, z), §(n + 1, z), ^{n + 2, z), . . . } and 
V: = Span{^'*(n, z), ^'*(n - 1, z), ^'*(n - 2, z), . . . }, 

with ^(n, z) and ^'*(n, z) the functions obtained by substituting for x in 
'I'(2;,z) and ^'*(x,z) respectively. As recalled above (see ()8.28|) . (|8.3U() and 
(jHinU), with W replaced by b{W)), the equations (lOHl) and (lO^ imply that 
^(n,z) and \^*(n,z) satisfy the bilinear identities i?(^(n, z), ^*(m, z)) = 
0, V n > m. Thus 

K* = V^, (8.35) 

where denotes the orthogonal of Vn with respect to B. 
Using H8.18() . we can write T{n,t) in 1)8. 10() as 

oo g 

T{n, t) = r(x, t) = exp(^ -^^iYl "^i) ^ 

i=i j=i 

det(x''^X'r^exp{t, Xj) - 6jf/-^exp{t, ^ij)) , (8.36) 

with X = g", and thus r(x, t) is a polynomial in x. Equations 1)8.3(1 and (|8.4() 
combined with the definition ()8.32() of ^'(x, z) and 'I'*(x, z), show then that 

K C z"r(z,t)-ip and V* C z'''T{qz,t)-^V, 
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which, because of (|8.H5|) . is equivalent to 

Clearly, Vn contains elements of orders n, n + 1, n + 2, . . . , and only those, 
meaning that the codimension of Vn in z'^t{z, t)~^V is equal to the degree (as 
a polynomial in z) of r(2;, t). Thus the flag V : • • • C Vn+i C Vn C Vn-i C • • • 
is a rational flag in the sense of Definition 18.21 and ^'(n, z) and ^*(n, z) are 
the (reduced) wave and adjoint wave functions associated with this flag. The 
lemma is established. □ 

Remark 8.10. Proposition 18 . 81 has been obtained by F.W. Nijhoff and O.A. 
Chalykh [22] in the special case when = q^i, by writing the (reduced) 
wave function as 

■^{x,z) = e'°'i°s'r ' detjld- (x ld + X)-^[X, Z]g-i{z Id + Z)"^}, 

with X and Z some g x g matrices, so that the g-commutator [X, Z]^-i is a 
rank 1 matrix. This argument is insufficient to deal with the (non-generic) 
situation that concerns us in Theorem 18. II They conjecture that all bispec- 
tral rank 1 commutative rings of difference operators can be parametrized 
in this way. In the limit q = 1, corresponding to rank 1 bispectral rings of 
differential operators, this was proved in a highly non-trivial paper by G. 
Wilson P- Iliev [2j has obtained independently of IH2] similar formulas. 
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